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Charging  of  a  Man  in  the  Wake  of  the  Shuttle 


I.  Introduction 

This  report  describes  &  portion  of  the  work  performed  by  S-CUBED  under 
the  contract  AFGL-F19628-82-C-0081  which  pertained  to  charging  of  an 
astronaut  performing  EVA.  The  report  covers  work  in  which  the  charging  of  the 
DMSP  satellites  is  studied  and  shown  to  be  understood  using  classical  models  of 
ion  collection.  The  physics  of  coupled  large  (shuttle)  and  small  (EMU)  object 
charging  in  strong  auroral  environments  which  was  developed  by  S-CUBED 
researchers  is  presented  here.  From  this  understanding,  the  POLAR,  NASCAP, 
and  MATCHG  codes  were  used  in  conjunction  to  calculate  the  charging  of  the 
EMU  in  the  auroral  environments  measured  by  the  DMSP  satellites  F6  and  F7. 
These  calculations  are  reported  here.  Much  of  the  material  presented  in  this 
report  has  been  presented  elsewhere  or  given  in  the  interim  report  presented  at 
the  Janurary  1986  Reno  AIAA  meeting. 

The  DMSP  satellite  has  been  observed  to  charge  to  potentials  as  large  as  679 
volts  negative  with  respect  to  the  local  plasma.  This  occurred  when  the  environ¬ 
ment  was  a  combination  of  an  intense  auroral  electron  precipitation  combined 
with  a  dropout  of  the  ambient  cold  plasma  density.  Previous  calculations  using 
the  POLAR  and  NASCAP  models  of  spacecraft  charging  demonstrated  that  the 
measured  environments,  when  used  as  input  to  the  codes,  predicted  charging  to 
the  levels  observed.  Calculations  of  charging  by  an  isolated  EMU  were  consistent 
with  the  DMSP  results.  Further  calculations  were  performed,  using  POLAR, 
which  showed  that  charging  by  shuttle  sized  objects  is  much  more  severe  than  for 
the  smaller  satellites.  An  environment  which  would  charge  DMSP  to  -172  Volts 
was  calculated  to  charge  'th^'  a  shuttle-sized  object  to  -3545  Volts.  However,  no 
theory  or  code  was  adequate  to  calculate  the  overall  and  differential  potentials  on 
the  EMU  during  EVA  in  the  vicinity  of  the  shuttle.  This  paper  discusses  such  a 
theory,  and  by  using  a  combination  of  computer  codes  is  able  to  present  calcula¬ 
tions  of  EMU  charging  during  shuttle  EVA. 

In  section  II.  an  analysis  of  the  charging  of  the  DMSP  satellite  is  presented 
using  the  MATCHG  code.  The  MATCHG  code  is  a  zero  dimensional  code  which 
has  an  orbit  limited  current  collection  model  and  the  sophistic ated^modefs  of  sur¬ 
face  secondary  and  backscatter  currents  contained  in  the  N  AS  CAP  code.  The 
MATCHG  analysis  of  DMSP  charging  matches  qualitatively  and  quantitatively 
the  measured  charging  without  any  free  parameters.  A  scaling  from  small  DMSP 

1 


rvaanvA  swjci 


KtcmMwasxvwa  urn tax  U9MAL  u*  kOMAUMfl 


size  objects  where  the  current  collection  is  orbit  limited  to  larger  shuttle  size 
objects  where  the  current  collection  is  space  charge  limited  is  presented.  Under 
DMSP  environments  the  shuttle  is  shown  to  charge  in  excess  of  3000  volts  nega¬ 
tive.  In  section  III.  the  physics  of  shuttle  EMU  charging  is  discussed.  Under  the 
effects  of  the  strong  auroral  environments  seen  by  the  DMSP  satellites,  the  shuttle 
is  shown  to  produce  a  high  energy  mono-energetic  ion  environment  for  the  EMU. 
The  energy  of  incoming  ions  is  approximately  the  shuttle  potential.  Differential 
charging  between  the  EMU  and  shuttle  and  between  different  materials  on  the 
EMU  is  shown  to  be  a  substantial  fraction  of  the  shuttle  potential.  In  section  IV. 
POLAR  calculations  of  the  shuttle  charging  in  the  DMSP  environments  are 
presented.  The  calculations  show  that  the  shuttle  charges  to  -3290  volts.  In  sec¬ 
tion  V.  MATCHG  calculations  of  EMU  materials  in  the  mono-energetic  environ¬ 
ment  determined  in  section  IV  are  presented.  These  calculations  show  that  charg¬ 
ing  of  the  EMU  to  over  6000  volts  negative  will  occur  and  differential  charging  of 
over  2000  volts  is  expected.  In  section  VI.  NAS  CAP  is  used  to  model  the  detailed 
charging  of  the  EMU  under  the  shuttle  induced  environments  determined  in  sec¬ 
tions  IV.  and  V..  The  NASCAP  calculations  confirm  the  MATCHG  calculations 
in  magnitude  of  charging  and  reveal  that  the  parts  of  the  suit  will  charge  much 
slower  than  the  rest  of  the  materials.  The  final  section,  VII,  gives  estimates  of 
the  energy  stored,  capacitively,  in  the  EMU  due  to  absolute  and  differential 
charging. 
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II.  Analysis  of  DMSP  Charging 

The  following  paper  titled  "ANALYSIS  OF  DMSP  CHARGING  AND 
IMPLICATIONS  FOR  POLAR  SHUTTLE  MISSIONS"  was  presented  at  the 
AIAA  shuttle  environment  interactions  workshop  held  in  Houston,  Tx  in 
November  1085.  The  paper  contains  an  analysis  of  the  observed  DMSP  F6  and 
F7  satellite  charging.  The  DMSP  satellite  charged  to  hundreds  of  volts  negative 
*  during  eclipse  conditions  of  high  auroral  electron  flux  and  low  ambient  ion  den¬ 

sity.  The  analysis  presented  shows  that  the  ion  current  collection  is  orbit  limited 
for  a  small  object  such  as  DMSP  and  will  become  space  charge  limited  for  shuttle 
site  objects.  The  MATCHG  code,  which  assumes  orbit  limited  current  collection 
and  has  NASCAP’s  models  of  surface  secondary  generation,  was  used  to  predict 
the  DMSP  charging  potentials  using  as  input  the  observed  environments  as 
presented  by  Gussenhoven  •t.oJ.I1!  The  MATCHG  predictions  of  the  DMSP 
charging  potentials  are  shown  to  agree  with  the  observed  charging  potentials  if 
the  ion  species  is  assumed  to  be  Oxygen.  For  a  hydrogen  ion  environment,  the 
predicted  charging  of  the  DMSP  satellite  is  shown  to  be  much  lower  than 
observed. 

I 
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ANALYSIS  OF  DMSP  CHARGING  AND  THE  IMPLICATIONS  FOR 
POLAR  SHUTTLE  MISSIONS 


I.  Katz,  M.  J.  Mandall,  C.  A.  Jongeward  J.  R.  Li  I  ley,  Jr. 

S- CUBED 

La  Jolla,  CA  92038 

WilliaatN.  Hall,  M.  S.  Gussenhoven 
Air  Fore*  Geophysics  Laboratory 
Hanscom  Air  Force  Base,  MA  01731 

The  high  levels  of  spacecraft  charging  seen  on  DMSP/F6  and  DMSP/F7 
resulted  from  the  combination  of  high  fluxes  of  auroral  electrons  and  a 
low  background  ambient  density.  Calculations  are  presented  using 
MATCHG,  a  material  charging  code,  which  show  the  observed  satellite 
charging  events  are  consistent  with  the  measured  auroral  environments. 
Full  three-dimensional  calculations  using  the  POLAR  code  with  these 
environments  show  that  a  shuttle  size  object  would  charge  to  thousands 
of  volts. 


Introduction 

High  levels  of  spacecraft  charging  have  been  observed  on  at  least 
two  satellites,  DMSP/F6  and  DMSP/F7,  in  840  km  polar  orbit. When 
charging  occurred,  the  satellites’  environments  were  characterized  by 
intense  fluxes  of  high  energy  auroral  fluxes  and  by  particularly  low 
densities  in  the  low  energy  ambient  plasma.  Calculations  using 
MATCHG^,  a  spacecraft  charging  model,  show  that  the  observed  levels  of 
charging  are  consistent  with  the  measured  environments.  The  model 
includes  secondary  and  backscattered  electrons,  which  msks  up  more  than 
half  of  the  positive  currents  to  the  satellite.  The  assumption  that  the 
ambient  plasma  consists  primarily  of  oxygen  ions  is  fundamental  to  the 
agreement  with  experiment.  Collection  of  ions  by  DMSP  is  shown  to  be 
orbit  limited. 
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For  larger  objects,  such  ss  th«  shuttle  or  an  sstronsut  near  the 
shuttle,  space  charge  effects  will  cause  substantially  greater  charging. 
Three-dimensional  calculations  using  the  POLAR  code  shoe  charging  of  the 
shuttle  in  excess  of  a  3000  volts  negative  will  occur  from  the  auroral 
events  seen  by  DMSP.  Even  at  higher  ion  dons i ties  more  common  to 
shuttle  altitudes,  charging  in  excess  of  1000  volts  negative  eill  occur. 

OMSP  Charging  Calculations 

The  OMSP  charging  calculations  were  performed  using  the  MATCHG 
code  for  the  three  environments  for  which  fits  to  the  high  energy 
auroral  electron  fluxes  were  provided  in  Reference  1.  As  a  worst  case 
scenario,  the  observation  with  the  most  intense  auroral  flux  was  chosen 
as  a  fourth  case  to  calculate.  The  auroral  flux  above  14  keV  was  0.4  x 
10®  (cm^  s-sr)“*  for  DMSP/F6  on  January  12,  1983.  In  the  absence  of  a 
spectra)  fit,  the  high  energy  electrons  were  assumed  to  have  a 
temperature  of  10.1  keV,  the  same  as  the  highest  intensity  case  with  a 
fit  reported.  For  this  case,  the  high  energy  auroral  electron  density 
used  was 


^observed 


(1) 


with  J0bservod/®  =  ®  4»10®/cm^/s-sr .  This  gives  the  value  for  n  of 
32/cm^  displayed  in  Table  2. 


The  MATCHG 
approximation  is 


code  assumes  orbit  limited  current  collection.  This 
valid  when  the  dominant  current  limiting  mechanism  is 
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the  conservation  of  angular  momentum  of  the  collactod  particles.  For  a 
DWSP  aiza  object,  tha  observed  charging  potantials  and  ambient 
anviromaant  produced  orbit  limited  conditions.  A  discuss* on  of  tha 
crossover  aiza  from  orbit  limited  to  space  charge  limited  collection 
will  be  deferred  to  tha  next  section.  In  this  section  we  will  assume 
current  collection  is  orbit  limited.  The  orbit  limiting  mechanism  is 
depicted  in  Figure  1  for  isotropic  and  flowing  plasmas.  For  the  flowing 
plasma  experienced  by  DMSP,  a  ram  ion  with  energy  ^ram  will  be  collected 
by  the  satellite  providing  angular  momentum  conservation, 

rsheath  vram  *  rsat  vsat 
and  energy  conservation, 

2 

v  . 

.  sat  . 

*^ram  "ion  2  *  ^sat 

can  be  satisfied  where 

,  _  1  mi on  y2 

^ram  2  q  ram 

Solving  these  simultaneous  equations  gives  the  radius  for  ion 
collection,  r,  the  current  to  the  satellite  is  then, 

s  2xr*  Jran, 

where  Jram  is  the  ram  ion  current  density.  Finally, 

=  2'rs,tJr.J1  *  f~) 
rram 

This  is  the  form  for  the  surface  current  used  by  the  MATCHG  code. 


Figure  1.  Orbit  limited  collection  from  a  flowing  plasma  with  velocity 
^ram- 


MATCHG  models  the  satellite  as  a  sphere  covered  with  a  single 
material  The  calculations  were  done  assuming  teflon,  since  teflon  beta 
cloth  is  a  typical  external  covering  of  low  altitude  spacecraft.  The 
external  surface  material  determines  the  secondary  and  backscattered 
electron  currents  which  cancel  over  three-fourths  of  the  incident 
auroral  electron  flux.  Table  1  shows  the  different  components  of  the 
current  before  charging  and  near  equilibrium. 
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The  top  half  of  Tablo  1  gives  tho  initial  charging  componanta 
corroaponding  to  tho  unchar god  DMSP  satellite  at  0.  volts.  Tho  lowor 
half  shows  tho  charging  components  for  tho  oqu *  I i br i urn  potential  reached 
by  tho  satellite.  Note  that  this  final  potential  represents  a  four 
order  of  Magnitude  reduction  in  the  charging  currents.  The  charging,  in 
this  case  to  *295  volts,  perturbs  the  electron  flux  only  slightly.  The 
ion  current  has  increased  sixty  fold  to  achieve  current  balance. 

Table  1.  MATCHG  Calculation  of  the  Charging  of  a  Teflon  Sphere 

Potential  0.  E+00  Volts 


Incident  Electron  Current  -1.05E-05 

Seconds  ry  E I ec trons  4 . 92E-06 

Back scattered  Electrons  2.98E-06 

Incident  Proton  Current  3.83E-08 

Secondary  Electrons  3.96E-10 

Bulk  Conductivity  Current  0.00E+00 


Net  Current  -2 .59E-06  amps/m^ 
Potential  -2.95«E02  Volts 


Incident  Electron  Current  -1.02E-05 

Secondary  Electrons  4.78E-06 

Backscattered  Electrons  2.90E-06 

Incident  Proton  Current  2.40E-06 

Secondary  Electrons  1.49E-07 

Bulk  Conductivity  Current  2.32E-10 


Net  Current  -6.36E-10  amps/m^ 

The  final  voltage  of  -302  corresponds  to  a  four  order  of  magnitude 
reduction  in  charging  current. 

The  charging  results  for  all  four  cases  are  shown  in  Table  3.  The 
environments  used  for  comparison  with  experiment  are  shown  in  Table  2. 
In  all  cases  the  ram  energy  was  taken  to  be  4.76  eV  corresponding  to 
oxygen  at  orbital  velocities. 
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Table  2.  Environments  Used  for  MATCHG  Calculations 


ni 

ne 

T.  (W) 

125 

3.9 

10.1 

190 

3.2 

14.4 

180 

4.9 

4.2 

537 

32.0 

10.1 

As  has  been  seen  in  studies  of  spacecraft  charging  at  geosynchronous 
orbit,  charging  requires  high  electron  temperatures  as  well  as  intense 
currents.  The  third  case  did  not  charge  because  a  4.2  keV  Maxwellian 
spectrum  generates  more  secondary  and  backscattared  electrons  on  a 
teflon  surface  than  incident  electrons. 

Table  3.  MATCHG  Charging  Calculations  for  Spacecraft  Potential 

Compared  with  Observation.  The  fit  to  the  F6  data  is 
discussed. 


Sate  1 1 i te 

UT 

0SC 

(Observed) 

^$C 

(Predicted) 

F7  11/26/83 

49843 

-317 

-295 

F7 

49844 

-140 

-284 

FT 

49845 

-  44 

-0 

F6 '  1/12/83 

35877 

-462 

-532 

*Data  fit  by  authors  (see  text) 

Further  studies  were  performed  to  evaluate  the  sensitivity  of  the 
charging  to  the  choice  of  spacecraft  external  material  and  to  the 
assumption  of  an  oxygen  plasma.  Table  4  shows  the  charging  levels  for 
four  different  spacecraft  materials.  Except  for  gold,  which  did  not 
charge,  the  three  other  materials  responded  similarly.  Because  gold  has 
a  high  atomic  number  it  is  relatively  easy  to  produce  secondary 
electrons  which  escape  preventing  charging.  Kapton,  which  has  been 
observed  on  scathaM  to  charge  strongly,  shows  charging  even  for  the 
mildest  environment.  The  level  of  charging  of  all  three  dielectric 
coatings  were  comparable  with  experiment. 
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Table  4.  Calculations  Using  MATCHG  for  Gold,  Kapton,  Teflon,  and  Indium 
Oxide  Spheres  Showing  the  Predicted  Charging  Levels  for  the 
Three  Spectral  Fits  from  Reference  1. 

Ram  Oxygen  Plasma 


Electrons 

Oxygen 

• 

N, 

T 

N. 

T 

[c»-3] 

[keV] 

[cm*3] 

[eV] 

Gold 

Kapt 

Tef  1 

Indo 

3.9 

10.1 

125 

4.76 

♦0.4 

-551 

-295 

-370 

3.2 

14.4 

140 

4.76 

-8.9 

-434 

-284 

-301 

4.9 

4.2 

180 

4.76 

♦1.4 

-163 

♦0.5 

-11.7 

The 

predicted  levels  of 

chargi ng 

change 

dramatically  if  the 

env i ronment  i s 

assumed  to 

be  dominated  by  hydrogen 

i  ons .  As 

seen  in  Eq. 

(3),  for 

sate  II i te  potent i a  1 s 

more  than 

a  few 

volts 

negative  the 

co 1 1 ec ted 

ion  < 

current  is  i 

inversely  proportional  to 

the  mass 

of  the  ion. 

To  collect  the  same  ion  current  in  a  hydrogen  plasma  requires  only  1/16 
the  potential  as  an  oxygen  plasma.  As  seen  in  Table  5,  the  equilibrium 
voltages  for  a  spacecraft  in  a  hydrogen  plasma  are  an  order  of  magnitude 
down  from  those  for  an  oxygen  plasma.  Even  if  the  ion  temperature  is 
assumed  to  be  1  eV,  the  predicted  potentials  for  all  materials  are  far 
below  those  observed  on  DMSP. 


Table  5.  Calculations  Using  MATCHG  Showing  the  Low  Level  Charging 
Predicted  for  a  Hydrogen  Background  Plasma  of  Hydrogen. 


Thermal  Hydrogen  Plasma 
Electrons  Hydrogen 


[cm*3] 

[keV] 

[cm-3] 

[eV] 

Gold 

Kapt 

Tef  1 

Indo 

3.9 

10.1 

125 

0.3 

♦0.4 

-35.8 

-17.8 

-22.7 

3.9 

10.1 

125 

1.0 

♦0.4 

-64.2 

-31.9 

-40.6 

3.2 

14.4 

190 

0.3 

-0.5 

-26.1 

-16.8 

-17.9 

3.2 

14.4 

190 

1.0 

-0.4 

-46.9 

-30.0 

-32.0 

4.9 

4.2 

180 

0.3 

♦1.4 

-  9.7 

♦  0.5 

-  0.6 

4.9 

4.2 

180 

1.0 

♦  1.4 

-17.2 

♦  0.5 

-  0.7 
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MATCHG  assumes  orbit  limited  collection  of  both  electron  and  ions. 
Orbit  limited  collection  of  the  energetic  electrons  is  easily  justified. 
For  the  ions,  orbit  limited  theory  will  also  be  valid  if  the  orbit 
limited  collection  radius  is  less  than  the  space  charge  limited  sheath 
distance  for  collection  of  the  same  ion  current  density  at  the  final 
voltage. 

The  orbit  limited  and  space  charge  limited  sheath  distances  for 
three  spherical  objects  with  radii  0.5  m,  1  m,  and  15  m  are  given  in 
Tabie  6.  The  first  two  radii  are  approximately  DMSP  size  objects  while 
the  15  m  radius  represents  a  shuttle  size  object.  The  orbit  limited 
sheath  radius  is  determined  by  angular  momentum  conservation 

rsheath  ”  Satellite 

A 

where  ^ram  is  the  ram  ion  kinetic  energy.  ^ram  »  m\n/2q.  For  oxygen  in 
an  800  km  orbit,  ^ram  is  4.76  volts.  Detailed  three-dimensional 
calculations  using  the  POLAR ^  code  indicate  that  the  ram  flux  produces 
an  approximate  isotropic  flux  with  magnitude  *TS3t^  Jram  =  4rrsai^  Jj . 
Using  this  value  for  the  flux  over  an  entire  spherical  sheath  boundary 
gives  the  space  charge  limited  radii  listed  in  Table  6.^ 


Table  6.  Comparison  of  the  Orbit  Limited  and  Space  Charge  Limited  Sheaths 
for  Two  Satellite  Size  Objects  and  One  Shuttle  Size  Object 


Jsc  W 

(Observed) 

Space  Charge 

n:  on 

Spacecraft  Orbit  Limited 

Limited  Sheath 

Sate) 1 

i  te  UT 

Rad i us 

Sheath  Radius 

Radius 

[cm”3] 

M 

Em] 

[m] 

F7 

49843 

-317 

125 

0.5 

4.2 

8 

1.0 

8.4 

11 

15.0 

125.0 

42 

F7 

49844 

-140 

190 

0.5 

2.8 

5  5 

1.0 

5.6 

7 

15.0 

85.0 

30 

F7 

49845 

-44 

180 

0.5 

1.6 

3.3 

1.0 

3.2 

4.6 

15.0 

48.0 

23 

F6 

35877 

-462 

537 

0.5 

5.1 

6.5 

1.0 

10.1 

8.6 

15.0 

152.0 
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Because  the  orbit  limited  sheath  is  within  the  range  of  space 
charge  limited  sheath  sizes,  an  orbit  limited  calculation  such  as  MATCHG 
is  appropriate  for  small  objects  such  as  DMSP.  However,  for  larger 
objects,  such  as  the  space  shuttle,  it  is  seen  that  the  current 
collection  is  strongly  space  charge  limited. 


>ace  Charge  Limited  Colloction 


As  is  seen  in  Table  6  and  Eq.  (3),  the  radius  for  orbit  limited 
current  collection  is  proportional  to  the  spacecraft  radius.  Hence, 
under  orbit  limited  conditions  a  spacecraft’s  equilibrium  potential  is 
independent  of  its  size,  because  the  surface  ion  flux  is  independent  of 
the  spacecraft  radius 
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2  2 

Surface  *  1  sheath'  r sat 


1  ♦ 


''sc'l 


ram 


This  result  is  valid  for  small  objects.  For  larger  objects,  the 
collected  current  produces  enough  space  charge  to  shield  the 
spacecraft’s  potential,  limiting  its  collection  ability.  This 
spacecharge  shielding  distance  introduces  a  new  length  scale  into  the 
physics.  So  the  size  invariance  of  the  equilibrium  potential  which 
existed  for  orbit  limited  collection  is  no  longer  true. 

For  spherical  objects  a  different  size  scaling  can  be  determined 
by  examining  Poisson's  equation  including  the  space  charge  contribution. 
For  completeness  this  scaling  argument  is  reviewed  in  Appendix  A.  The 
result,  Eq.  (A. 3)  in  the  appendix,  is  2.  Spacecraft  with  sizes  Lj,  Lj 
will  have  equilibrium  potentials  related  by  =  (Lj/I^)*^.  This 
effect  is  reflected  in  Table  6  where  the  space  charge  limited  sheath 
distance  scales  much  slower  than  the  spacecraft  radius.  For  oxample,  in 
the  auroral  environments  seen  by  DMSP,  a  two  meter  radius  object  (a  2 
meter  size  object  is  just  beginning  to  be  affected  by  space  charge 
shielding)  charged  to  a  few  hundred  volts.  In  contrast,  a  shuttle  size 
object  (approximately  15  meters  in  radius)  would  charge  according  to  Eq. 
(A. 3)  to 


shutt I e 


=  (XS/2)4/3  =  15 


2m 


(4) 


A  shuttle  size  object  will  charge  than  an  order  of  magnitude  more  than  a 
2  meter  size  object  in  the  same  environment. 


Shuttle  Col  lection  and  Astronaut  Char 


The  calculations  oV  auroral  charging  on  shuttle  size  objects  were 
accomplished  using  the  POLAR  computer  code.  The  POLAR  code  is  a  three- 
dimensional  code  which  computes  the  potentials  on  and  about  a  large 
object  (with  respect  to  the  Debye  length)  accounting  for  space  charge, 
spacecraft  wakes,  secondary  electrons  and  backscatter.  The  auroral 
environments  used  are  given  in  Table  7.  The  first  environment 
corresponds  to  the  observation  on  DMSP/F7  at  49843  UT.  The  second 
environment  is  the  "strongest*  auroral  environment  which  was  used 
previously  with  the  MATCHG  calculations,  except  the  density  of  the  cold 
ions  has  been  increased  to  104/cm®,  a  more  typical  value  for  a  300  km 
polar  orbiting  shuttle. 

Table  7.  Charging  of  a  Shuttle-size  Object 

n.(c»*3)  kTa(I.V)  n j (cm-3)  #SC(V)  Ws».l|) 

3.9  10.1  125  -3290  (-172) 

32  10.1  104  -1306  (-23) 

The  model  used  in  the  POLAR  calculations  had  the  approximate  size 
and  shape  of  the  actual  shuttle  and  is  shown  in  Figure  2.  Teflon  was 
chosen  as  a  covering  material  for  the  model  because  it  has  similar 
secondary  characteristics  as  the  shuttle  insulating  tiles.  The 
computational  space  was  a  cubical  mesh  of  size  3.5  meters  and  an  overall 
grid  size  of  (120  meters).^]  Although  this  model  is  a  crude  rendition 
of  the  actual  shuttle,  the  bulk  charging  results  are  reliable  because 
the  size  attributes  and  material  properties  are  comparable. 


The  results  of  the  POLAR  cs I  cut  at  ions  along  with  the  environment! 
are  given  in  Table  7.  The  column  ^sc  is  the  POLAR  result  for  the 
equi I ibrium  potential  of  the  shuttle  size  object.  The  final  column, 
given  for  comparison,  is  the  potential  to  which  a  QMSP  size  object  would 
have  charged  in  the  same  environment.  The  charging  of  shuttle  size 
objects  is  seen  to  be,  as  expected,  much  more  severe  than  for  smaller 
objects. 

Figure  3  shows  the  potential  contours  from  the  POLAR  calculation. 
The  view  is  from  the  side  with  the  shuttle  nose  at  the  top  of  the 
figure.  The  motion  of  the  shuttle  is  to  the  left.  Because  of  the 
high  charging  potential,  the  sheath  is  spherical  and  far  from  the 
shuttle.  The  flowing  ions  from  the  ram  direction  were  able  to  reach  the 
wake  side  of  the  shuttle.  This  produced  a  spherical  distribution  of 
ions  masking  much  of  the  ram  wake  effects.  For  this  auroral  environment 
the  shuttle  charged  more  than  3000  volts  negative. 

The  potential  contours  for  the  second  POLAR  calculation  are 
presented  in  Figure  4.  For  this  higher  cold  ion  density  the  shuttle 
charging  was  less  severe  but  still  greater  than  1500  volts  negative. 
The  sheath  is  much  closer  to  the  shuttle  and  larger  ram  wake 
differential  charging  is  predicted,  although  the  difference  was  less 
than  250  V  since  this  simple  model  assumed  the  shuttle  was  covered  with 
a  single  material.  Experience  from  spacecraft  charging  at 
geosynchronous  orbit  shows  differential  potentials  arise  primarily  from 
the  responses  of  different  materials. 

It  is  important  to  realize  that  an  astronaut  performing  EVA  during 
this  auroral  event  would  charge  to  similarly  high  potentials  depending 
on  his  location  with  respect  to  the  shuttle.  Far  from  the  shuttle, 
being  a  small  object,  an  astronaut  would  charge  to  potentials  similar  to 
DMSP.  Near  the  shuttle,  however,  charging  to  high  levels  comparable  to 
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the  shuttle  would  occur.  For  tho  low  background  density  case,  an 
astronaut  would  experience  little  ram  wake  differences.  Ram  wake 
differences  increase  with  increasing  background  dens  ity. 
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gure  2.  Shuttle  mode  used  in  the  POLAR  calculations. 
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The  POLAR  calculations  presented  showed  that  shuttle  size  objects 
will  charge  to  thousands  of  volts  in  the  DMSP  observed  environments. 
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Potential  Scali 


Taking  a  plasma  with  temperature  T,  and  density  N,  for  Poisson’s 
aquation 


-na 


If  rsheath 
Ur#  r 

4 


2 


(A.l) 


Equation  (A.l)  is  valid  for  potentials  much  larger  than  the  plasma 
temperature  because  the  repelled  species  has  been  ignored.  Let 
^0*sat*l)  denote  a  solution  to  Eq.  (A.l)  with  boundary  conditions 
appropriate  to  a  satellite  of  radius  R. 


#(Ksat4  *  *sat 
t Rshaath*  ® 

It  is  straightforward  to  check  by  substitution  into  Eq.  (A.l)  that 
#(r.b)  .  b4/3*(r/b,  1) 

is  a  solution  to  Poisson's  equation  corresponding  to  a  satellite  b  times 
as  large,  with  a  sheath  b  times  farther  away,  that  is, 


1  T 

^rsheath 

Ur#(r,b) 

r 

.  4 

f(bR5tt,b)  =  b4/VRslt.i) 

"0R,h..th^  *  b4/3"(R.h..th'1)  *  0 


(A.  2) 
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Henee  as  the  object  size  increases  the  object  potential  l(bRsat»b),  must 
increase  to  keep  the  surface  flux  constant.  By  Eq.  (A. 2),  given  two 
spacecraft  with  sizes  Lj  und  L2  their  equilibrium  potentials  wi 1 1  be 
related  by 


(A.3) 
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III.  The  Physics  of  Shuttle  EMU  Charging 


Background 

Figures  1  and  2  show  the  POLAR  models  of  the  shuttle  and  EMU  used  dur¬ 
ing  the  calculations.  The  combined,  multiple  grid  model  is  shown  in  Figures  3.  It 
is  immediately  apparent  that  the  EMU  is  extremely  small  compared  with  the 
shuttle.  Indeed  examination  of  the  space  potentials  for  the  combined  system,  Fig¬ 
ure  4,  shows  that  the  electric  fields  due  to  the  EMU  are  not  felt  a  few  EMU  radii 
from  the  EMU.  The  theory  presented  below  shows  that  ion  trajectories  are 
modified  only  slightly  by  the  presence  of  the  EMU.  It  is  also  shown  that  for  the 
EMU  to  collect  the  ion  current  density  needed  to  balance  auroral  fluxes,  the  EMU 
charges  to  much  higher  potentials  when  near  the  shuttle  than  when  isolated.  The 
quantatative  theory  predicts  differential  charging  potentials  comparable  with 
overall  charging  values. 


Surface  Charging  Physics 

The  primary  cause  for  differential  charging  on  an  object  as  small  as  the  EMU 
is  the  differing  secondary  and  backscatter  currents  of  the  various  surface  materi¬ 
als.  The  charging  current  to  each  surface  is  the  difference  between  the  incident 
electron  current  and  the  outgoing  secondary  and  backscattered  currents.  For  a 
given  distribution  of  electron  energies  there  can  be  a  large  range  of  net  charging 
currents  depending  upon  the  surface  material.  For  the  intense  energetic  auroral 
electron  fluxes,  some  materials  such  as  gold  emit  about  one  electron  for  each 
incident  electron;  others,  such  as  kapton,  emit  only  a  third  the  incident  electrons. 
The  surfaces  charge  negatively  until  the  difference  between  incident  and  emitted 
electrons  is  made  up  by  ions.  Thus  for  objects  with  varying  surface  materials  the 
ion  current  density  must  also  vary  along  the  surface.  This  requires  the  surface 
potentials  to  vary  accordingly.  The  question  of  how  the  surface  current  density 
varies  as  a  function  of  voltage  distribution  has  been  addressed  previously  for 
current  collection  by  high  voltage  solar  arrays.  That  work  is  extended  here  to  cal¬ 
culate  ion  collection  by  a  very  small  object  in  the  presence  of  a  large  ion  collecting 
object. 
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Figure  1.  POLAR  Model  of  Shuttle. 
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Figure  2.  POLAR  and  NASCAP  Model  of  EMU. 


Figure  3.  Multi-grid  Combined  POLAR  Model  of  EMU  and  Shuttle. 


28 


The  currents  which  produce  surfacing  charging  are  composed  of  four  parts, 


where  and  are  the  auroral  electron  and  collected  plasma  ion  current 

densities  incident  on  the  surface,  JHt  is  the  secondary  electron  current  emitted  due 
to  ion  and  electron  impact  and  JM  is  the  electron  backscatter  current.  The  last 
three  currents  each  contribute  a  positive  current  to  the  surface.  Ion  induced  secon¬ 
dary  electrons  are  ignored  in  these  calculations  because  of  a  lack  of  experimental 
data  on  the  materials  of  the  EMU  and  Shuttle.  Photo-currents  have  not  been 
included, hence  these  calculations  can  be  considered  as  eclipse  calculations. 

NASCAP  and  POLAR  use  a  secondary  emission  model  which  is  based  on  the 
deposition  of  energy  in  the  material  by  the  incident  electron  and  the  ability  of  the 
electrons  liberated  from  atoms  by  this  energy  to  reach  the  surface  of  the  material. 
The  six  parameters  needed  to  model  the  secondary  electron  emission  are  given  in 
the  equations  below.  A  detailed  description  of  this  model  is  given  in  the  NASCAP 
users  manual^2!. 


Where, 


R=blE0ni+b2E2ni 


The  parameters  n1  ,  bx,  nt  ,  are  entered  as  material  properties  in  the  object 
definition  of  NASCAP  and  POLAR.  For  convenience,  the  parameters  C  and  A 
are  computed  by  the  code  from  the  inputed  values  of  the  maximum  yield  value 
and  energy  value.  The  values  used  for  the  EMU  were  determined  by  fitting  ab- 
initio  calculations  of  secondary  electron  emission  by  Ashley  et.al.  !3!  to  the 
NASCAP  form  given  in  the  above  equation.  These  values  are  given  below  by 
material. 
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LEXAN 

# 

PROPERTY 

INPUT  VALUE 

1 

DIELECTRIC  CONSTANT 

3.50E+00  (NONE) 

3.50E+00  (NONE) 

2 

THICKNESS 

1.00E-02  METERS 

1.00E-02  MESH 

3 

CONDUCTIVITY 

1.00E-16  MHO/M 

1.00E-16  MHO/M 

4 

ATOMIC  NUMBER 

5.00E+00  (NONE) 

5.00E+00  (NONE) 

5 

DELTA  MAX>COEFF 

1.80E+00  (NONE) 

2.36E+01  (NONE) 

6 

E-MAX  >DEPTH**-1 

2.00E-01  KEV 

5.59E-02  ANG-01 

7 

RANGE 

7.15E+01  ANG. 

4.29E+01  ANG. 

8 

EXPONENT>  RANGE 

6.00E-01  (NONE) 

5.13E+02  ANG. 

9 

RANGE>  EXPONENT 

2.90E+02  ANG. 

6.00E-01  (NONE) 

10 

EXPONENT 

1.77E+00  (NONE) 

1.77E+00  (NONE) 

11 

YIELD  FOR  1KEV  PROTONS 

4.55E-01  (NONE) 

4.55E-01  (NONE) 

12 

MAX  DE/DX  FOR  PROTONS 

1.40E+02  KEV 

1.40E+02  KEV 

13 

PHOTOCURRENT 

2.00E-05  A/M* *2 

2.00E-05  A/M* *2 

14 

SURFACE  RESISTIVITY 

1.00E+16  OHMS 

8.85E+04  V-S/Q 

WHITE  PAINT 

# 

PROPERTY 

INPUT  VALUE 

1 

DIELECTRIC  CONSTANT 

3.50E+00  (NONE) 

3.50E+00  (NONE) 

2 

THICKNESS 

1.00E-02  METERS 

1.00E-02  MESH 

3 

CONDUCTIVITY 

5.90E-14  MHO/M 

5.90E-14  MHO/M 

4 

ATOMIC  NUMBER 

5.00E+00  (NONE) 

5.00E+00  (NONE) 

5 

DELTA  MAX>COEFF 

2.10E+00  (NONE) 

4.06E+01  (NONE) 

6 

E-MAX  >DEPTH**-1 

1.50E-01  KEV 

8.74E-02  ANG-01 

7 

RANGE 

7.15E+01  ANG. 

4.29E+01  ANG. 

8 

EXPONENT>  RANGE 

6.00E-01  (NONE) 

5.52E+02  ANG. 

9 

RANGE>  EXPONENT 

3.12E+02  ANG. 

6.00E-01  (NONE) 

10 

EXPONENT 

1.77E+00  (NONE) 

1.77E+00  (NONE) 

11 

YIELD  FOR  1KEV  PROTONS 

4.55E-01  (NONE) 

4.55E-01  (NONE) 

12 

MAX  DE/DX  FOR  PROTONS 

1.40E+02  KEV 

1.40E+02  KEV 

13 

PHOTOCURRENT 

2.00E-05  A/M* *2 

2.00E-05  A/M*  *2 

14 

SURFACE  RESISTIVITY 

1.00E+13  OHMS 

8.85E+01  V-S/Q 

Table  1.  EMU  material  parameters  used  in  the  NASCAP  calculations. 


TEFLON 


# 

PROPERTY 

INPUT  VALUE 

1 

DIELECTRIC  CONSTANT 

2.00E+00  (NONE) 

2.00E+00  (NONE) 

2 

THICKNESS 

1.00E-02  METERS 

1.00E-02  MESH 

3 

CONDUCTIVITY 

1.00E-16  MHO/M 

1.00E-16  MHO/M 

4 

ATOMIC  NUMBER 

7.00E+00  (NONE) 

7.00E+00  (NONE) 

5 

DELTA  MAX>COEFF 

2.40E+00  (NONE) 

3.00E+01  (NONE) 

6 

E-MAX  >DEPTH**-1 

3.00E-01  KEV 

1.46E-01  ANG-01 

7 

RANGE 

4.54E+01  ANG. 

1.81E+01  ANG. 

8 

EXPONENT>  RANGE 

4.00E-01  (NONE) 

1.24E+02  ANG. 

9 

RANGE>  EXPONENT 

7.00E+01  ANG. 

4.00E-01  (NONE) 

10 

EXPONENT 

1.77E+00  (NONE) 

1.77E+00  (NONE) 

11 

YIELD  FOR  1KEV  PROTONS 

4.55E-01  (NONE) 

4.55E-01  (NONE) 

12 

MAX  DE/DX  FOR  PROTONS 

1.40E+Q2  KEV 

1.40E+02  KEV 

13 

PHOTOCURRENT 

2.00E-05  A/M* *2 

2.00E-05  A/M* *2 

14 

SURFACE  RESISTIVITY 

1.00E+16  OHMS 

8.85E+04  V-S/Q 

ALUMINUM 

# 

PROPERTY 

INPUT  VALUE 

1 

DIELECTRIC  CONSTANT 

1.00E+00  (NONE) 

1.00E+00  (NONE) 

2 

THICKNESS 

1.00E-03  METERS 

1.00E-03  MESH 

3 

CONDUCTIVITY 

-1.00E+00  MHO/M 

-1.00E+00  MHO/M 

4 

ATOMIC  NUMBER 

1.30E+01  (NONE) 

1.30E+01  (NONE) 

5 

DELTA  MAX>COEFF 

9.70E-01  (NONE) 

9.18E+00  (NONE) 

6 

E-MAX  >DEPTH**-1 

3.00E-01  KEV 

3.00E-02  ANG-01 

7 

RANGE 

1.54E+02  ANG. 

1.23E+02  ANG. 

8 

EXPONENT>  RANGE 

8.00E-01  (NONE) 

3.87E+02  ANG. 

9 

RANGE>  EXPONENT 

2.20E+02  ANG. 

8.00E-01  (NONE) 

10 

EXPONENT 

1.76E+00  (NONE) 

1.76E+00  (NONE) 

11 

YIELD  FOR  1KEV  PROTONS 

2.44E-01  (NONE) 

2.44E-01  (NONE) 

12 

MAX  DE/DX  FOR  PROTONS 

2.30E+02  KEV 

2.30E+02  KEV 

13 

PHOTOCURRENT 

4.00E-05  A/M**2 

4.00E-05  A/M**2 

14 

SURFACE  RESISTIVITY 

-1.00E+00  OHMS 

-8.85E-12  V-S/Q 

Table  1.  Continued 
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KAPT 

# 

PROPERTY 

INPUT  VALUE 

1 

DIELECTRIC  CONSTANT 

3.50E-f00  (NONE) 

2 

THICKNESS 

1.00E-02  METERS 

3 

CONDUCTIVITY 

1.00E-16  MHO/M 

4 

ATOMIC  NUMBER 

5.Q0E+00  (NONE) 

5 

DELTA  MAX>COEFF 

1.80E+00  (NONE) 

6 

E-MAX  >DEPTH**-1 

2.00E-01  KEV 

7 

RANGE 

7.15E+01  ANG. 

8 

EXPONENT>  RANGE 

6.00E-01  (NONE) 

9 

RANGE>  EXPONENT 

2.90E+02  ANG. 

1° 

EXPONENT 

1.77E+00  (NONE) 

11 

YIELD  FOR  1KEV  PROTONS 

4.55E-01  (NONE) 

12 

MAX  DE/DX  FOR  PROTONS 

1.40E+02  KEV 

13 

PHOTOCURRENT 

2.00E-05  A/M**2 

14 

SURFACE  RESISTIVITY 

1.00E+16  OHMS 

Table  1.  Continued 


3.50E+00  (NONE) 
1.00E-02  MESH 
1.00E-16  MHO/M 
5.00E-S-00  (NONE) 
2.36E+01  (NONE) 
5.59E-02  ANG-01 
4.29E+01  ANG. 
5.13E+02  ANG. 
6.00E-01  (NONE) 
1.77E+00  (NONE) 
4.55E'01  (NONE) 
1.40E+02  KEV 
2.00E-05  A/M* *2 
8.85E+04  V-S/Q 
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The  graph  below  gives  the  secondary  electron  emission  curves  for  the  five 
materials  of  the  EMU.  Note  that  Kapton  and  Lexan  have  the  same  secondary 
electron  characteristics.  This  is  because  the  secondary  properties  of  Lexan  were 
not  known.  Kapton  was  chosen  as  a  model  of  Lexan  because  they  are  both  insu¬ 
lating  plastics. 
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Figure  5. 
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Secondary  Electron  Emission  Coefficients  for  the  EMU  Materials 
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Backscattered  electrons  are  those  emitted  from  the  surface  with  energies 
above  50  eV.  Their  energy  distribution  is  usually  peaked  close  to  the  primary 
incident  energy  and  they  may  be  considered  as  reflected  electrons.  POLAR  and 
NASCAP  employ  an  empirical  representation  for  the  backscattered  flux  depend* 
ing  only  on  the  atomic  number  of  the  material  and  the  incident  angle.  This  form 
is  suggested  by  experimental  data  end  Monte  Carlo  data^,  and  is 

17 

where 

-E 

e*-[/o,(S0£)e(r-.0S)^|i-+8(£-l)l|.l«  * 

where  E  is  in  KeV,  Z  is  the  atomic  number  of  the  material,  and  9  is  a  step  func¬ 
tion. 


The  material  properties  used  to  model  the  shuttle  charging  are  given  below. 


WHITE  SPUTTERED  TILES 

# 

PROPERTY 

INPUT  VALUE 

1 

DIELECTRIC  CONSTANT 

2.00E+00  (NONE) 

2.00E+0Q  (NONE) 

2 

THICKNESS 

1.00E-02  METERS 

1.00E-02  MESH 

3 

CONDUCTIVITY 

1.00E-16  MHO/M 

1.00E-16  MHO/M 

4 

ATOMIC  NUMBER 

7.00E+00  (NONE) 

7.00E+00  (NONE) 

5 

DELTA  MAX>COEFF 

1.60E+00  (NONE) 

5.91E+00  (NONE) 

6 

E-MAX  >DEPTH**-1 

6.00E-01  KEV 

1.14E-02  ANG-91 

7 

RANGE 

4.54E+01  ANG. 

1.81E+01  ANG. 

8 

EXPONENT>  RANGE 

4.00E-01  (NONE) 

3.24E+02  ANG. 

9 

RANGE>  EXPONENT 

2.00E+02  ANG. 

4.00E-01  (NONE) 

10 

EXPONENT 

1.62E+00  (NONE) 

1.62E+00  (NONE) 

11 

YIELD  FOR  1KEV  PROTONS 

4.55E-01  (NONE) 

4.55E-01  (NONE) 

12 

MAX  DE/DX  FOR  PROTONS 

1.40E+02  KEV 

1.40E+02  KEV 

13 

PHOTOCURRENT 

2.00E-05  A/M* *2 

2.00E-05  A/M* *2 

14 

SURFACE  RESISTIVITY 

1.00E+16  OHMS 

8.85E+04  V-S/Q 

Table  2.  Material  properties  for  White  Sputtered  tiles. 
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Parameter*  5  through  10  in  the  table  above  were  chosen  to  fit  the  measured 
secondary  emission  coefficients  for  the  shuttle  white  tiles  as  measured  by  Yang 
Data  for  sputtersd  tiles  was  taken  because  shuttle  tiles  are  most  certainly 
sputtered  by  Oxygen  bombardment  in  the  orbit.  The  secondary  electron  emission 
curve  for  this  material  is  given  in  figure  6.  The  triangles  show  the  experimental 
values  of  Yang  «(.•/. 
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Figure  6.  Secondary  Electron  Emission  Coefficient  for  the  White  Sputtered 

Tiles 
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An  area  which  needs  further  investigation  is  the  effect  of  secondary  electrons 
generated  by  ion  impact.  Although  NASCAP  and  POLAR  account  for  this  posi¬ 
tive  surface  current,  the  models  in  the  code  rely  on  experimental  data  which  is  not 
accurately  known.  A  preliminary  examination  of  recent  data  suggests  that  ion 
generated  secondaries  can  be  a  substantial  portion  of  the  positive  current  to  the 
surface  for  a  shuttle  size  object  in  the  environments  seen  by  the  DMSP  satellites 
F6  and  F7.  This  is  based  on  rough  estimates  for  the  ion  produce  secondaries  for 
aluminum  taken  from  a  recent  publication  by  Alonso  et.a/J6)  and  should  be  con¬ 
sidered  tentative.  Increased  ion  generated  secondaries  would  decrease  the  shuttle 
charging  potentials  over  the  values  calculated  in  this  analysis,  however  shuttle 
charging  of  at  least  a  thousand  volts  is  still  expected  and  differential  charging 
ratios  predicted  will  still  be  qualitatively  correct. 
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Differential  Voltages  and  Charging 

In  this  section  we  will  obtain  an  estimate  of  the  differential  charging  induced 
on  an  EMU  performing  EVA  in  the  wake  of  the  orbiter.  The  physical  situation  to 
be  im'estig ated  is  therefore  a  very  small  object  (an  EMU  of  about  .5m  in  radius) 
near  a  very  large  object  (an  orbiter  of  about  15m  in  radius).  This  is  shown  in  the 
POLAR  code  plot,  figure  2.  As  will  be  shown,  the  electric  fields  a  few  meters 
from  the  EMU  are  determined  by  the  orbiter  potential.  The  ion  collection  situa¬ 
tion  for  the  EMU  is  thus  similar  to  a  probe  in  a  monoenergetic  ion  beam  where 
the  beam  energy  is  the  spacecraft  potential.  If  the  EMU’s  material  properties  are 
similar  to  the  nearby  orbiter  then  the  EMU  will  charge  to  the  shuttle  potential 
and  will  collect  an  ion  flux  equal  to  that  of  the  orbiter.  However,  if  the  EMU’s 
material  surface  secondary  and  backscatter  properties  are  different,  then  to  reach 
current  equilibrium  the  EMU  must  charge  to  a  potential  different  from  the  shut¬ 
tle.  This  will  increase  or  decrease  the  collected  ion  flux.  As  will  be  shown,  a  frac¬ 
tional  change  in  differential  potential  between  the  orbiter  and  the  EMU  is  needed 
to  cause  a  corresponding  fractional  change  in  the  collected  current.  Because 
differences  in  secondary  properties  between  different  materials  can  cause  net  col¬ 
lected  current  to  vary  by  50  to  75  percent,  differential  charging  of  the  EMU  in  the 
wake  of  the  orbiter  can  be  expected  to  be  on  the  same  scale  as  the  average  orbiter 
charging. 


Local  Electric  Fields 

In  this  section  we  will  estimate  the  magnitude  of  the  electric  field  perturba¬ 
tions  caused  by  a  small  object  in  the  vicinity  of  a  much  larger  charged  object. 
The  situation  is  depicted  in  Figure  7.  The  large  sphere  represents  the  orbiter  and 
the  small  sphere  represents  the  EMU. 

Let  L  be  the  radius  of  the  large  sphere,  1  the  radius  of  the  small  sphere  and  d 
their  separation.  The  potential  of  the  large  sphere  is  $  and  the  smaller  spvere 
differs  from  this  potential  by  A$.  The  region  we  are  investigating  is  L  >>  d  >>  1 
and  The  crossover  point,  r,  where  an  ion  becomes  influenced  by  the 

differential  potential  can  be  estimated  as  the  point  where  the  electric  field  from 
the  differential  potential  equals  the  electric  field  from  the  large  sphere.  This  elec¬ 
tric  field  is  due  to  the  potential  of  the  charged  large  sphere  (orbiter)  and  is  given 
by. 


El-9 


L 

(r+i)J 


where  r  u  measured  from  the  surface  of  the  large  sphere.  Similarly  the 
differential  electric  field  due  to  the  small  sphere  being  at  a  different  potential  than 
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the  larger  sphere  can  be  estimated  as 


l 

(r-d)* 


As  shown  in  Figure  7,  r  is  the  distance  from  the  surface  of  the  larger  sphere  and 
taken  on  a  radius  which  passes  through  the  centers  of  both  spheres.  This  gives 
an  upper  bound  on  the  strength  of  the  field  due  to  the  differential  charging 
because  the  dipole  nature  of  the  fields  is  being  ignored  and  its  inclusion  would 
result  in  a  faster  drop  in  the  differential  voltage  field. 

This  gives  for  the  crossover  distance  in  the  limit  L  >>d  >>  I, 


r-*\/Z7 

Using  values  for  L  and  1  corresponding  to  an  orbiter  and  an  EMU,  i.e., 

L—15mcttn 


we  obtain  that  for  the  extreme  case  where  the  differential  charging  equals  the 
potential  of  the  orbiter,  that  at  distances  greater  than  3  meters  from  the 

EMU,  the  collected  ions  will  not  be  influenced  by  the  EMU.  At  this  distance 
from  the  EMU  and  the  orbiter  the  space  potential  is  approximately  the  orbiter 
potential,  ie. 

since, 


Differential  Current  Collection 

We  have  shown  that  a  EMU  charged  to  a  potential  A*  different  from  a 
nearby  large  orbiter  disturbs  the  local  electric  fields  only  a  few  meters  from  the 
EMU.  Because,  as  shown  above,  the  space  potential  at  the  crossover  region  is 
approximately  the  orbiter  potential,  we  can  estimate  the  current  collected  by  the 
EMU  using  the  following; 

l)The  presence  of  the  large  shuttle  creates  an  environment  for  the  EMU  consist¬ 
ing  of  a  mono-energetic  beam  of  ions  with  energy  equal  to  the  orbiter  potential. 
This  is  because  the  ions  are  accelerated  to  the  orbiter  potential  before  they  are 
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substantially  effected  by  the  EMU  electric  fields,  as  discussed  above. 

2)Ion  collection  by  the  small  EMU  in  the  wake  will  be  orbit  limited,  with  a 
current  density  for  surfaces  at  the  orbiter  potential  the  same  as  the  average 
orbiter  ion  current  density.  This  is  because  a  small  probe  in  a  ion  beam  collects 
according  to  orbit  limited  theory. 

The  ion  collection  is  then  given  by, 

Hence,  once  we  have  determined  the  charging  of  the  orbiter  to  get  *  we  can  in  a 
second  calculation,  using  the  formulation  discussed  above,  determine  the  ion 
environment  at  the  EMU,  and  from  this  the  EMU  charging.  This  last  step  can  be 
approximated  using  MATCHG  to  estimate  AJ,  the  variation  in  charging  currents 
due  to  surface  secondary  and  backscatter  differences,  which  using  the  above  equa¬ 
tion  then  would  give  the  differential  voltages,  A$,  needed  to  produce  these  varia¬ 
tions.  A  more  accurate  NASCAP  calculation  can  be  performed  using  the  effective 
ion  environment  the  shuttle  would  produce.  These  two  approaches  will  be 
presented  in  sections  V  and  VI. 


43 


IV.  POLAR  Code  Calculations 


The  POLAR  code  was  used  to  determine  the  coarse  charging  that  an  orbiter 
EMU  pair  would  produce  in  an  auroral  charging  environment.  The  model  used 
for  the  calculations  is  shown  if  figure  8.  Although  this  is  a  crude  model  of  the 
orbiter,  the  over  all  size  and  shape  is  similar.  The  material  covering  the  model 
has  secondary  properties  equivalent  to  the  shuttle  white  tile  material  and  a  small 
patch  on  the  wake  side  is  made  of  teflon  to  simulate  the  EMU. 

The  auroral  environment  used  in  the  POLAR  calculations  was  taken  from 
the  uts  to  the  observed  DMSP  environments  given  by  Gussenhoven  et.  al.  This 
was  10.1  KeV,  for  the  auroral  electron  temperature  and  3.9/em3  for  the  density. 
The  cold  ion  density  was  taken  to  be  the  observed  density  of  125/cm3. 

POLAR  code  calculations  were  performed  for  a  range  of  voltages  from  -1000 
volts  to  -4500  volts.  From  these  calculations  a  current  verses  voltage  relationship 
was  determined.  A  plot  of  I  Vs  V  for  the  shuttle  shown  in  Figure  9.  Also  shown 
on  this  figure  is  the  net  electron  current  to  the  shuttle  surface.  This  current  is 
composed  of  three  terms 

^ntt m  ^incident +  ^tecondary^haekacaUtr 

The  last  two  terms  in  the  above  equation  contribute  a  positive  current.  The  point 
of  intersection  of  the  two  curves  is  the  steady  state  charging  potential  of  the  shut¬ 
tle,  For  reasons  discussed  in  section  III,  ion  generated  secondaries  have  been 
neglected  in  this  calculation.  From  Figure  9  it  is  clear  that  shuttle  charging  is 
reduced  as  ion  generated  secondary  emission  is  increased.  Figure  10  shows  con¬ 
tour  plots  for  the  charged  shuttle  at  the  -3290  volt  equilibrium  potential. 
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Figure  10.  POLAR  Potential  Calculation  of  the  Charged  Shuttle 


V.  MATCI  r  Calculations  of  EMU  Charging 

The  environment  used  in  the  MATCHG  calculations  was  identical  to  that  for 
the  POLAR  calculations  and  discussed  in  section  IV.  The  hot  electron  spectrum 
was  a  single  Maxwellian  with  temperature  10.1  KeV  and  a  density  of  3.9/cm3. 
The  ion  temperature  and  density  are  chosen  so  that  the  slope  of  the  ion  current 
collection  verses  voltage  is  the  correct  orbit  limited  value  derived  in  section  III. 
This  turned  out  to  be  .5 2/cm3.  The  temperature  is  chosen  .0048  KeV  correspond¬ 
ing  to  the  ram  energy  of  ambient  oxygen  ions.  With  this  choice  of  plasma  tem¬ 
perature  and  density  MATCHG  will  correctly  model  the  orbit  limited  charging  of 
the  materials. 

Table  2  shows  the  charging  potentials  and  charging  currents  for  8  different 
materials  typical  of  the  shuttle  and  EMU.  As  is  seen,  variations  of  3000  volts 
over  the  EMU  and  between  the  EMU  and  the  shuttle  are  expected  due  to 
differences  in  material  secondary  and  backscatter  properties.  A  summary  of  the 
material  charging  is  shown  below  for  the  8  materials. 


Material 

Net  Electron  Current 

Potential 

Kapton 

2.86  uA/M**2 

-6290  volts 

White  Tiles 

1.15 

-3290  * 

(sputtered) 

Teflon 

1.56 

-4090 

Aquadag 

2.67 

-6000  * 

White  Non 

4.04 

-2700 

Conducting  Paint 

Gold 

0. 

.4  * 

Alum 

2.20 

-6800 

*Not  a  suit  material. 

Table  2.  Charging  potentials  for  selected  materials.  The  column 
labeled  Net  Charging  Current  is  sum  of  the  incident  electron 
current,  secondaries  and  backscatter. 
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VI.  NASCAP  Calculations 

In  this  section  we  will  present  two  calculations  of  the  differential  charging 
expected  on  the  EMU  in  the  wake  of  a  Polar  orbiting  shuttle.  NASCAP  calcula¬ 
tions  were  performed  using  the  EMU  model  shown  in  Figure  2.  The  material  pro¬ 
perties,  such  sc  secondary  and  backscatter  emissions,  were  chosen  to  closely  model 
those  of  the  actual  EMU.  Figures  11a  through  Ilf  show  the  EMU  material  con¬ 
struction  as  modeled  in  the  POLAR  and  NASCAP  calculations.  As  in  the 
MATCHG  calculations,  the  hot  electron  temperature  and  density  where  take 
from  fits  to  DMSP  data  given  by  Gussenhoven.  The  cold  ion  environment  was 
chosen  similar  to  the  MATCHG  calculations  where  the  temperature  is  taken  to  be 
low  to  simulate  the  mono-energetic  nature  of  the  charging  environment  of  the 
EMU,  and  the  density  is  chosen  to  give  the  proper  orbit  limited  current  verses 
voltage  characteristic. 

The  results  of  the  NASCAP  calculations  are  given  in  Table  3.  Reading  down 
the  columns  shows  the  convergence  of  the  NASCAP  charging  calculation  indicat¬ 
ing  that  the  final  few  iterations  are  well  converged  for  all  but  the  ground  conduc¬ 
tor.  Fox  the  aluminum  components  and  the  suit  ground,  the  differential  charging 
time  scale  is  tens  of  seconds  which  was  longer  than  the  NASCAP  calculations 
were  run.  Figure  13  shows  the  potentials  of  the  suit  ground  (ALUMINUM)  and 
the  kapton  gloves  as  a  function  of  time.  The  aluminum  is  still  charging  while  the 
kapton  has  reached  the  current  balance  predicted  by  the  MATCHG  calculations. 
This  long  charging  time  scale  results  from  the  relatively  small  current  collection 
area  for  the  suit  ground.  The  charging  time  scale  for  the  suit  ground  can  be 
estimated  using  the  incident  electron  current  to  the  aluminum  patches  and  the 
capacitance  between  the  suit  surface  and  the  underlying  ground. 
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where  Cf-<0  ***  is  the  suit  capacitance,  Jtlte  is  the  net  electron  current  to  the 

surface,  and  is  the  exposed  area  of  suit  ground  (aluminum)  and 

A(H£-5.8m2  is  the  suit  surface  area.  From  the  MATCHG  calculations  this  gives 

Jtl  e—  2.3>*-4r  *  charging  rate  of  in  agreement  with  the  -16-*^- 

in*  of  see  see 

from  Table  3.  Over  100  seconds  are  required  to  reach  the  -6800  volts  expected 

from  the  MATCHG  calculations  The  MATCHG  calculations  showed  that  the 

aluminum  surfaces  should  charge  to  more  than  6000  volts  negative  given  sufficient 
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time.  Figures  14  and  15  show  potential  contours  about  the  EMU  for  this  charg¬ 
ing  calculation.  As  seen,  potentials  on  the  face  plate  of  the  astronaut  are 
thousands  of  volts  higher  than  on  the  suit  material  of  the  EMU. 

Differences  between  MATCHG  and  NASCAP  result  from  NASCAPs  proper 
inclusion  of  three  dimensional  effects.  Current  balance  is  effected  by  leakage 
currents  to  underlying  conductors  and  surface  conductivity.  In  addition,  adjacent 
surfaces  which  are  charging  differently  due  to  secondary  property  differences,  can 
effect  each  other's  charging  characteristics  becuase  of  electric  field  effects  on  secon¬ 
dary  emission. 
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POTENTIAL  IN  KILO-VOLTS 

TIMS 

S 

C 

D 

8 

r 

Q 

H 

4.1843 

-.444 

-.994 

-.993 

-.493 

-.493 

-.493 

-.441 

1.71 -OS 

•149 

•149 

•149 

•149 

•1.34 

•149 

•144 

1. 41-03 

-3.04 

•3.04 

•3.03 

-3.04 

•3.04 

•3.04 

•3.04 

3.1842 

-ITS 

•3.73 

•3.73 

-3.73 

•173 

•3.73 

-3.47 

3.1843 

-341 

441 

441 

449 

-3.39 

*3.31 

•344 

1.3841 

•443 

443 

443 

4.13 

4.14 

443 

•4.04 

348-01 

-4.41 

4.41 

4.41 

4.37 

4.39 

4.41 

•4.31 

348-01 

•440 

440 

440 

440 

4.31 

4.30 

•4.34 

4.38-01 

-4.38 

44S 

443 

4.30 

-4.33 

444 

-4.34 

(.38 -01 

•4.40 

4.90 

4.00 

4.39 

4.33 

4.41 

-4.33 

6.38-01 

•4.43 

4.43 

•4.93 

4.34 

443 

4.43 

•443 

T.38-01 

•4.70 

4.70 

•4.70 

444 

443 

4.70 

•4.33 

(48-01 

-4.78 

4.73 

4.74 

447 

443 

4.73 

-443 

948-01 

•4.79 

4.79 

4.79 

447 

4.33 

4.40 

-4.33 

1.08+00 

•4.44 

4.34 

443 

4.34 

4.33 

4.44 

•443 

148+00 

-4.94 

4.94 

4.93 

443 

4.34 

4.93 

•441 

148+00 

4.04 

4.04 

4.03 

443 

4.33 

-3.04 

-441 

1.38+00 

-3.13 

4.13 

4.13 

4.34 

4.33 

•313 

•4.30 

3.08+00 

-S.33 

443 

440 

-443 

4.34 

•3.33 

•4.20 

348+00 

-440 

440 

4.34 

443 

444 

-3.30 

-440 

348+00 

•347 

447 

443 

443 

447 

-347 

-4.19 

3.48+00 

-5.44 

4.44 

4.43 

443 

444 

-3.44 

-4.19 

3.0B+00 

-341 

440 

4.44 

441 

4.34 

•341 

-4.19 

348+00 

-344 

444 

444 

440 

449 

-347 

-4.1( 

348+00 

-4.0 

4.43 

449 

4.30 

449 

-3.43 

-4.18 

3.38+00 

-3.47 

4.47 

4.44 

4.19 

4.40 

•3.47 

-4.1( 

4.08+00 

-4.71 

4.71 

449 

4.19 

4.40 

-3.72 

-4.17 

448+00 

-4.74 

4.74 

4.73 

4.14 

4.40 

•4.74 

-4.17 

4.38+00 

4.90 

4.40 

4.77 

4.14 

4.41 

-3.(0 

-4.17 

4.IB+00 

4.44 

4.44 

4.30 

4.14 

4.41 

-3.(4 

-4.18 

(.08+00 

4.47 

447 

4.44 

4.17 

4.41 

-3.(7 

-4.18 

(48+00 

4.90 

4.90 

4.47 

4.17 

4.43 

•3.90 

•4.18 

S4E+00 

4.93 

4.93 

4.90 

4.17 

4.43 

•3.93 

•4.18 

(48+00 

4.94 

4.94 

4.93 

4.14 

4.43 

-3.94 

•4.15 

e.oa+00 

4.99 

4.99 

4.95 

4.14 

4.43 

-3.99 

-4.13 

448+00 

4.01 

4.01 

4.97 

4.14 

4.43 

•4.01 

•4.13 

4.(8 +00 

4.03 

4.03 

4.01 

4.15 

4.44 

-4.03 

•4.15 

748+00 

4.09 

4.09 

4.05 

4.13 

4.44 

-4.09 

•4.14 

7.48+00 

4.13 

4.13 

4.03 

-448 

4.43 

•4.13 

•4.14 

448+00 

4.15 

-8.15 

4.11 

4.14 

4.48 

•4.15 

-4.14 

448+00 

-8.17 

4.17 

4.13 

4.17 

4.48 

-4.17 

-4.14 

Table  3  Results  of  NASCAP  calculation  of  EMU  charging  in  the  near  shuttle 
environment.  The  letters  B,C,D,E,F,G  refer  to  locations  on  the  EMU  shown  in 
Figure  12 
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SlMPACt  CCll  MATCOIAL  COMPOSITION  AS  «Itt»  MOM  TNC  POSIWC  T  OtMCTION 

roo  t  vacucs  octwcn  i  am  it 

MATfOIAC  LCQCNO 


Figure  11  (a)  NASCAP  Material  Plots  of  the  EMU. 
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1.  S.  K.  7.  9.  11.  13.  IK.  17.  l».  21.  23.  2K  .  27.  29.  31.  33. 

-2- 


Figure  11.  (b)  NASCAP  Material  Plots  of  the  EMU 


53 


MATERIAL  LECENO 


1 

Lai* 


SURFACE  CELL  MATERIAL  COMFOSITIOM  AS  VIEWED  FROM  THE  POSITIVE  I  DIRECTION 
FOR  I  VALUES  OETREEN  t  AMO  }1 


>•  ~  ■  *•  *■  *•  T-  »•  t.  If.  li.  11.  1J.  i«.  ii,  11.  it. 


Figure  il.  (c)  NASCAP  Material  Plots  of  the  EMU. 
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■  ■ 


i.  a.  i.  4.  i.  «.  7.  a.  a.  it.  u.  12.  13.  14.  it.  la.  i  ? . 


Figure  SI.  (d)  NA5CAP  Material  Plots  of  the  EMU. 
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Potential  (Kilovolts) 


Figure  13. 


Time  History  of  Charging  of  the  Suit  Ground  and  Kapton  Gloves 


)}.!• 


I 


i 

s 

I 

s 

31, M 
2MI 

27. M 

i 

26.00 

23.60 

21.00 

10.00 

17.00 

16.00 

13.00 

S 

11.00 

9.00 

7.00 

6.00 

3.00 

1.00 . 

12346070  9  10  11  12  13  14  16  16  17 

X-AXIS 

-514E*03  <  CONTOUR  LEVELS  <  -1.43E*93  AT  POTENTIAL  INCREMENTS  OF  OF  »  2.00E«O2 

1.00  <X<  17.00,  1.00  <Z<  33.00,  CUTPLAHE  OFFSET  V*  9.00 


Figure  14  Potential  Contours  Around  Charge  EMU. 
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I 


-e.l7E*»3  <  CONTOUR  LEVELS 

!.§•  <V<  37.M, 


Y-AXIS 

•  X , 4  3E«03  AT  POTENTIAL  INCREMENTS  OF  DP  *  6 

l.fR  <Z<  33.00,  CUTPLANE  OFFSET  X«  9.00 


Figure  15.  Potential  Contours  Around  Charge  EMU. 
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VD.  Stored  Electrical  Energy  in  the  EMU 

Energy  is  stored  in  the  EMU  through  capacitance  to  infinity  and  to  underly¬ 
ing  conductors.  The  quantity  of  energy  in  each  of  these  modes  depends  on  the 
voltages  and  capacitances  as, 

where  6  V  is  the  voltage  difference  between  the  suit  surface  material  and  the  suit 
ground,  where  V  is  the  surface  voltage  relative  to  infinity,  and  C is  the  capaci¬ 
tance  to  infinity.  These  capacitances  can  be  estimated  by  approximating  the 
EMU  as  a  conducting  spherical  ball  with  radius  R  and  covered  with  teflon  of 
thickness  d.  The  capacitance  to  infinity  can  then  be  determined  by  the  definition 
of  capacitance. 


Q~CV 

and  the  relationship  between  charge  and  potential  on  a  sphere, 

This  gives  for  the  capacitance  of  a  sphere, 


C00»4«0R 

-102  pF  (1  meter  object) 

The  capacitance  to  the  underlying  conductor  (suit  ground)  Cf  is  much  larger  and 
can  be  estimated  assuming  a  parallel  plate  capacitor  of  area  the  surface  area  of 
the  EMU.  The  capacitance  of  a  parallel  plate  capicator  is, 


C, 


where  k  is  the  dielectric  constant  which  we  will  call  1  for  these  estimates,  A  is  the 
surface  area,  and  d  the  thickness  of  the  insulating  suit  material  .  For  a  spherical 
EMU  the  surface  area  is  A=4nR2  so  the  capacitance  to  suit  ground  is 

n  a  r2 
C,  =  47Tf0— 


which  for  a  1  meter  radius  EMU  with  1  centimeter  of  insulating  material  results 
in 


Cs=ltfpF 


Hence  Cr=-^-  Coc=100C(X, 


so  for  comparable  voltages  much  more  energy  is  stored 
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The  MATCHG  calculations  showed  that  the  dominant  material,  teflon, 
would  charge  to  approximately  >4000  volts  and  that  the  face  plate  material  would 
charge  to  *6000  volts,  and  the  suit  ground  to  about  -6800  volts.  Hence, 
4000*  and  f  V--1000*.  Using  these  values  EM  and  Eey  can  be  estimated  as, 

E'-.l  Joule, 

These  calculations  are  upper  bounds  on  the  energy  involved  in  a  discharge 
because  they  assume  all  the  stored  charge  is  transferred  in  the  discharge.  In  real¬ 
istic  discharges  only  a  fraction  of  the  charge  (depending  on  the  geometry  and  the 
type  of  discharge)  contributes  to  the  discharge. 

These  estimates  can  be  improved  using  the  results  of  the  NAS  CAP  calcula¬ 
tion.  The  NASCAP  code  computes  the  capacitance  to  infinity  for  a  three  dimen¬ 
sional  object  as  part  of  its  charging  algorithm.  For  the  EMU  NASCAP  gives  the 
effective  radius  of  flOSii/  thus  C00-fi8p/\  and  £oo-.3i0~3Jo*i««. 

Two  other  important  points  should  be  made  at  this  point.  First,  grounding 
the  astronaut  would  increase  the  arcing  hazard  by  fixing  the  suit  ground  at  the 
shuttle  ground.  This  would  increase  SV  by  a  factor  of  three  and  hence  the  stored 
energy  by  ten  fold.  Secondly,  The  charging  of  the  suit  ground  happens  on  time 
scale  which  is  much  slower  than  the  other  charging  times  because  of  the  relatively 
small  amount  of  exposed  ground  conductor.  Hence,  for  early  times  (as  seen  in  the 
NASCAP  charging  history  shown  section  VI)  the  differential  voltage,  SV  is  over 
twice  the  steady  state  value  of  800-1000  volts. 


VII.  Stored  Electrical  Energy  in  the  EMU 

Energy  U  stored  In  the  EMU  through  capacitance  to  infinity  and  to  underly- 
ing  conductors.  The  quantity  of  energy  in  each  of  these  modes  depends  on  t^e 
voltages  and  capacitances  as, 

where  SV  is  the  voltage  difference  between  the  suit  surface  material  and  the  suit 
ground,  where  V  is  the  surface  voltage  relative  to  infinity,  and  cw  is  the  capaci¬ 
tance  to  infinity.  These  capacitances  can  be  estimated  hy  approximating  the 
EMU  as  a  conducting  spherical  ball  with  radius  R  and  covered  with  teflon  of 
thickness  d.  The  capacitance  to  infinity  can  then  be  determined  by  the  definition 
of  capacitance, 

QmCV 

and  the  relationship  between  charge  and  potential  on  a  sphere, 


4«0  R 


This  gives  for  the  capacitance  of  a  sphere, 


Coomiwt0R 


-10 *pF  (imtterobjeet) 

The  capacitance  to  the  underlying  conductor  (suit  ground)  Cf  is  much  larger  and 
can  be  estimated  assuming  a  parallel  plate  capacitor  of  area  the  surface  area  of 
the  EMU.  The  capacitance  of  a  parallel  plate  capicator  is, 


c»"*‘o-7 


where  k  is  the  dielectric  constant  which  we  will  call  1  for  these  estimates,  A  is  the 
surface  area,  and  d  the  thickness  of  the  insulating  suit  material  .  F*r  a  spherical 
EMU  the  surface  area  is  A-iwR2  so  the  capacitance  to  suit  ground  is 


C -4*cr 


which  for  a  1  meter  radius  EMU  with  1  centimeter  of  insulating  material  results 


Ct-lO*PF 

R 

Hence  Cf*-jCoe-lOOCoe,  so  for  comparable  voltages  much  more  energy  is  stored 
in  the  capacitance  to  suit  ground  than  in  the  capacitance  to  infinity. 
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The  MATCHG  calculation*  showed  that  the  dominant  material,  teflon, 
would  charge  to  approximately  *4000  volts  and  that  the  face  plate  material  would 
charge  to  *0000  volts,  and  the  suit  ground  to  about  *0800  volts.  Hence, 
4000#  and  4K—1000*.  Using  these  values  Em  and  Etv  can  be  estimated  as, 

Egm  .lJoulat 

These  calculations  are  upper  bounds  on  the  energy  involved  in  a  discharge 
because  they  assume  all  the  stored  charge  is  transferred  in  the  discharge.  In  real¬ 
istic  discharges  only  a  fraction  of  the  charge  (depending  on  the  geometry  and  the 
type  of  discharge)  contributes  to  the  discharge. 

These  estimates  can  be  improved  using  the  results  of  the  NASCAP  calcula¬ 
tion.  The  NASCAP  code  computes  the  capacitance  to  infinity  for  a  three  dimen- 
B:onal  object  as  part  of  its  charging  algorithm.  For  the  EMU  NASCAP  gives  the 
effective  radius  of  R,ffm  .608 A/  thus  and 

Two  other  important  points  should  be  made  at  this  point.  First,  grounding 
the  astronaut  would  increase  the  arcing  hazard  by  fixing  the  suit  ground  at  the 
shuttle  ground.  This  would  increase  SV  by  a  factor  of  three  and  hence  the  stored 
energy  by  ten  fold.  Secondly,  The  charging  of  the  suit  ground  happens  on  time 
scale  which  is  much  slower  than  the  other  charging  times  because  of  the  relatively 
small  amount  of  exposed  ground  conductor.  Hence,  for  early  times  (as  seen  in  the 
NASCAP  charging  history  shown  section  VI)  the  differential  voltage,  6V  is  over 
twice  the  steady  state  value  of  800-1000  volts. 
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